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A field-theoretical representation is suggested for the electron global density of states distribution 
function V{y) in extended disordered conductors. This opens a way to study the complete statistics 
of fluctuations. The approach is based on a functional integration over bi-local functions \P ri ,r 2 
instead of the integration over local functions in the usual functional representation for moments 
of physical quantities. The formalism allows one to perforin the disorder averaging and to derive 
an analog of the usual nonlinear a- model - a "slow" functional of a supermatrix field Q rijI . 2 (r) ~ 
^r,T! ° ^r 2 ,r- As an application of the formalism, the long-tail asymptotics of V[y) is derived. 
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A powerful theoretical tool for studying localization 
phenomena in disordered media is based on the func- 
tional integral representation of the Green's functions. 
This representation is the starting point for the deriva- 
tion of an effective slow action - the nonlinear cr-model 
0,0- The renormahzation group (RG) analysis of this 
model describes the scaling of the averaged conductance 
< g > with the system size L in accordance with the one- 
parameter scaling hypothesis [3J. Since the discovery of 
mesoscopic fluctuations (see review Q) the importance 
of studying distribution functions has been realized. For 
systems of dimension d < 1, distribution functions of 
some physical parameters can be found with the use of a 
"zero-dimensional" supersymmetric cr-model or transfer- 
matrix formalism (see references in the books |2|,|J]). For 
extended systems (d > f), distribution functions of some 
local quantities - like cigcnfunction amplitudes or the lo- 
cal density of states (DOS) - have been obtained |5| by 
means of the saddle-point approach |6( to the supersym- 
metric cr-model. However, no suitable field-theoretical 
representation for the global DOS and conductance distri- 
bution functions, V(y) and V(g), in systems of dimension 
2 < d < 4 has been known up to now; the information 
about these functions has been extracted from calcula- 
tion of all the moments < v n > and < g n > pj . A field 
theory for distribution functions is certainly in demand. 

In the present paper we develop such a theory for V(y) 
and derive a slow effective functional Eqs.J3J)-(Q} that de- 
termines the characteristic function T v {s) : 

V u (s) =< exp [sv/v] >= J V[Q] exp [F + F s ] , (1) 

where < . . . > denotes the averaging over the disorder 
potential U; v =< v >, and the global DOS, v, is deter- 
mined for a particular realization of the disorder by 

v{E) = i/(2irL d )Tr{G R (E) - G A (E)} . (2) 

Here G R( - A \E) = [E - H ±id}- 1 are exact retarded and 
advanced electron Green's functions for the Hamiltonian 
H = Hq + U; the trace Tr is taken over space arguments. 



For the unitary ensemble (spinless electrons in the 
presence of a time-inverse symmetry breaking magnetic 
field) the "free" action F is given by 

F ° = T J dT Str {^( V ^( r )) 2 - 2 <5A z Q(r)} ; (3) 

D = vpl/d is the diffusion coefficient; I >> pp 1 is the 
mean free path, pp and vp are the Fermi momentum and 
velocity; S is a dephasing rate. The "source" action F s is 

F s = -s J drStr{0<">(r)A,*Q(r)} , (4) 

where S = stA/(2it), t = l/v F , and A = \/(vL d ). 

In contrast to the usual cr-models, supermatrices Q in 
Eqs.0J,J3l, an d 0) depend on three spacial arguments: 
Q(r) = {Qnr 2 ( r )}; the dependence on r is smooth. The 
symbol "Str" denotes a combination of the usual super- 
trace operation str{Af} = tr{M bb } - tr{M ff } (b and / 
refer to boson and fermion components) and the sum- 
mation over the (lower) spacial arguments: Str{M} = 
s t r {-^riri }• The summation runs over sites of an 
auxiliary lattice of spacing a; in the continuum limit 
(a — > 0) the summation is replaced for J dri/a d . 

In Eqs.lPJ-Q, Q is a 4 x 4 matrix in the kinetic R-A 
(retarded - advanced) space and in the b-f superspace; 
Aj (i = x,y,z) are the Pauli matrices acting in the R- 
A space; k = diag(l,— 1) is a supermatrix in the b-f 
superspace. The supermatrix Q obeys the constraint 

Q 2 (r)=I, (5) 

where / is the unity in both the intrinsic matrix space and 
the coordinate space. The operator (r) is defined as 

OM 2 (r)=a 2ri 5[r-(r x +r 2 )/2]; (6) 

C = [1 + 4t 2 £ 2 ]~ 1 , where £ is the operator of kinetic 
energy counted from the Fermi energy Ep. 

The slow functionals Eqs.© and 0} allow one to study 
the complete statistics of fluctuations v in disordered con- 
ductors. Earlier, the way of studying V(g) and V(y) was 
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to calculate all the moments < g n > and < v n >. This 
formidable problem has been attacked in with the use 
of an extended cr-model. The matrix variable Q l,J (r) of 
this model has additional indices (i, j = 1, . . . , n) to ac- 
count for coupling (at the averaging over the disorder) of 
primary fields ^(r) and ^(r) used in the standard func- 
tional integral representation of each term (i = 1 , . . . , n) 
of a product g n or v n . The action of the extended a- 
model includes all the terms of the expansion in powers 
of IVQ and ojtQ, while the usual cr-model corresponds to 
the "hydrodynamical" approximation, i.e. to the lowest 
non-zero terms ~ (VQ) 2 and ~ u>Q; ui is the frequency 
of an external field. The formal reason of importance of 
short distances has been revealed by the RG analy- 
sis of the extended cr-model: vertices with 2n gradients 
have a positive anomalous dimension ~ (n 2 — n). This 
anomaly leads to an "explosion" of high cumulant mo- 
ments, so that the distribution functions Vig) and V(y) 
possess logarithmically normal (LN) tails [7j : 

V{y) oc exp [-(l/4u) In 2 [Sv/utA]] , (7) 

where 5v = v — v and u = In (oo /a) (oo is a bare value 
of the conductivity a for the orthogonal ensemble). The 
same mechanism leads also to the long-time asymptotics 
of the averaged conductance < g(t) > 0; it has been 
interpreted in terms of the LN distribution of the current 
relaxation times t$: 

ViU) oc exp [-(l/4u) In 2 (t 4 /r)] . (8) 

Assuming that the asymptotics Eq.JJJ) is caused by a 
single "quasi-localized" state of a small energy width 
~ l/ttf, with the corresponding contribution 5v ~ t^ 
to the glob al DOS, one may obtain Eq.0 directly from 
Eq.lJSJ |lCtlll|. On the other hand, no similar reasoning 
connecting Eq.JSJ) with the LN asymptotics of the con- 
ductance distribution function has been proposed yet. 

The distributions Eq.JJJ) and Eq.JBJ have been ob- 
tained 00 with the use of the sophisticated RG analysis 
restricted formally to d = 2 + e case. Later on, an al- 
ternative approach to study the asymptotics of < g(t) > 
has been suggested in 0. Based on the saddle-point 
analysis of the usual (non- extended) supersymmetric cr- 
model for the averaged conductance, this approach re- 
produces Eq.0) at d=2 and gives results for other di- 
mensions |g, |l2j . The problem of short distances in this 
approach was treated by introducing a cutoff r* ~ I [lflj . 

It is tempting to confirm in a similar way the RG re- 
sult Eq.JJJ) for the global DOS (as well as for the con- 
ductance) distribution function. However, this task has 
remained unfeasible because of the absence of a work- 
ing field-theoretical representation for V(y) (and V(g)). 
Below we shall achieve this goal with the use of the de- 
veloped theory, Eqs.© and (@J. 

For completeness, we briefly comment on a seemingly 
simple but hardly useful representation for V(y). It is 



based on the integration in Eq.Q over supermatrices 
Q l > 3 (r) which have only one spacial argument but an in- 
finite number of colors: 1 < i,j < AT; AT — * oc. The 
free action Fo is still given by Eq.® (with Str — > str), 
and F s = -(a/W) J str{A z kQ(r)}dr / L d . However, the 
advantage of the supersymmetry - a finite rank of the 
Q-matrix - is lost and it is not clear how to treat non- 
perturbative regimes; e.g. there is no apparent way to 
find solutions to the saddle-point equation V(QVQ) = 
[sA/(2nAfD)][A z k,Q(r)] and the corresponding action, 
which would behave properly in the limit AT — ► co. 

The formalism leading to the effective field theory of 
interest, Eqs.® and is based on the following rep- 
resentation of V v (s) in the form of a functional integral 
over bi-local primary superfields $ rirj and f rir2 : 

V V ( S ) = J £>[# , *] < e *[*A4*-S(Z3*+*£S)] > (9) 

Here the superfield ^ ri r2 is defined as a column of 
complex (S R ( A )) and Grassmann (£- R (- A )) fields ^ = 
(S R ,S A ;£ R ,£ A y and * = *+L, where the matrix L = 
diag(Lb,L/) = diag(A z ,7); the transposition "t" does 
not touch space variables, while the Hermitian conjuga- 
tion is defined with their interchange. An explicit form 
of the scalar product is ^S/M'i 1 = ^r 2 r 1 A4 rir2;r3T4 1 i T4T3 
(with the summation over repeated indices), where 

M ri r 2 ;r 3 r 4 = S T2 r 3 T[E -H + iSA z /2] rir4 . (10) 

In Eq.©, *6rir 2 = BS rir2 S a b (and similarly for B) has 
only the boson component; the constant vectors B and 
B in the R-A space are specified by their direct product 
BoB = il = A z + iA y . The choice of the matrix Vl is 
dictated by the requirements of the factorizability of 
and convergence of integrals in the slow functional (de- 
tails will be presented in an extended publication). We 
choose B = diag(l, 1) and B = diag(l, —1)'. Integrating 
over \&,^I> in Eq.lJHJ, we find: 

V v {s) =< exp[-Strln{M} - isStr{Af _1 nni,II}] > (11) 

Here n& = (1 + k)/2 is the projector on the boson sector 
of the b-f superspace; n riI . 2;r3r4 = S Tir2 S r3r4 . The first 
term in the exponent vanishes due to the supersymme- 
try (matrix A4 is diagonal in the b-f superspace) while 
the second one equals —sv/v, thus providing the desired 
equality T , v {s) =< exp[— sv/v] >. 

The primary representation Eq.® allows one to 
average the term exp [— i T , ^.rK-^r/] over the 
Gaussian disorder potential U and to use the 
Hubbard-Stratonovich (HS) decoupling of the term 
(^ , rir \E , rri )(\E , r2 r v I' rr2 ) with an auxiliary matrix field 
Q ri r 2 (r) of the usual symmetry = LQL (see |2|). Af- 
ter this transformation, the bilinear over 'J and ^ part of 
the action Eq.Q will be modified in the following way: 
M -> M(Q) = Mo+iQ/2, where M is given by Eg . lfTU|) 
with H = H ; Q rir2;r3r4 = S rirA Q T2r3 (ri). Integration 
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over ^ leads to Ea. ljTT|l where M is replaced by M(Q) 
and the averaging goes over the Gaussian action for the 
HS field Q. A (primary) saddle-point approach to the 
modified Eq. reveals a manifold Eq. © where the ac- 
tion has a degenerate extremum perturbed by the <5A 2 /2 
term and by the "source" term (oc s). Restricting to 
only the "transverse" variations of Q (obeying Eq.©) 
and keeping only the lowest powers of VQ (this hydro- 
dynamical approximation corresponds to a non-extended 
cr-model), we obtain Eq.© for V„(s) with the "free" ac- 
tion Fq Eq.© and the "source" action F s given by 

F s = -s J drStr{OM(r)(A»+»A 9 )Q(r)(l+*)}.(12) 

With the use of the parametrization Q: Q — A Z (I + 
iP)(I - iPy 1 (with A Z P = -PA A, we find that the 
"Az-part" of the action Ea. (|12|l is even in P, while the 
"A y -part" is odd. The diagrams generated by the A y -part 
of Ea. l|12|l are similar to those which would be caused by 
higher terms of the gradient expansion of Strln{. . .}, if 
we allowed for effects of the extended cr-model. In the 
considered hydrodynamic approximation all such terms 
should be ignored. Also, having in mind the assumed 
smallness of the source term, we shall keep in Eq. l(T2^) 
only the part (oc k) which breaks down the symmetry 
between the boson and fermion components of the fr- 
action (if not this part, the integral Eq.© would be just 
unity for any s). Omitting in Ea. (|12|l the unity term of 
the projector (1 + k) and the "A y -part", we arrive at the 
announced result, Eq.©. 

To check the agreement with the usual diagrammatic 
technique, we treat F s as a perturbation and calcu- 
late ln[7 3 „(s)] as a series in powers of s. This calcula- 
tion can be performed with the use of the mentioned 
parametrization of Q in terms of the non-constrained field 
P. The leading contribution to the term of a given or- 
der in s is made by diagrams with a minimal number of 
loops formed by lines associated with propagators ( "dif- 
fusons") < PP > of the field P. For instance, for the 
the first and second order terms of ln[V^(s)} we find — s 
and s 2 A 2 /(27r) 2 ^ q _^ [_Dq 2 ] -2 , respectively. From here 
we obtain the average < vjv >— 1 and the variance 
< (5v/v) 2 >= A 2 /{2T: 2 )J2q=ioi D( l 2 ^ 2 coinciding with 
the results of the diagrammatic technique jl4| . The lat- 
ter expression for the variance < (5v/v) 2 > is determined 
by one-loop diagrams. Higher cumulant moments of &v 
are determined by diagrams with higher number of loops 
0. Using Eqs.Q and ©, we have checked a non-trivial 
"accidental" cancellation of two-loop contributions to the 
third-order cumulant < (5v/D) 3 > c , so the latter is de- 
termined by three-loop diagrams, in accordance with 0. 

Most interesting applications of the developed formal- 
ism, Eqs.© and ©, are expected in the non-perturbative 
regime. In the present paper we implement the idea of the 
approach and consider non- trivial "secondary" saddle- 
point configurations of the action F = Fq + F s . Looking 



for an extremum of F with respect to variations of the 
field Q (obeying Eq.©): Q -> Q + SQ, SQ = [Q, e], we 
obtain the following "secondary" saddle-point equation: 



V[Q(r)VQ(r) 



[A,fcOM(r),Q(r)]. (13) 



Below we restrict the analysis to a particular sub- 
manifold £q of the saddle-point manifold £: matrices 
Qrir 2 ( r ) £ £o are characterized by a smooth depen- 
dence on the "center-of-mass" coordinate (ri+r2)/2 (i.e., 
\QIV'Q\ -C 1) and by a small difference |ri — r2 1 ~ I. 

It is convenient to re- write Eqs.© and © m the 
Wigner representation over the spacial indices: A TlT2 — 
(« rf /i d )E P ^(r';p)exp[ip(r 1 -r 2 )]; r' = (n + r 2 )/2. 
Omitting for brevity the <5-term in Eq.©, we have: 



Fn = 



^£str{^[V r Q(r;r';p)] 2 }, (14) 



/dr 
-^^C p str{A z fcQ(r;r;p)}, (15) 
p 

where C p = [1 + 4t 2 ( 2 ]~ 1 ; ( p = vp{p — pp). In the con- 
sidered hydrodynamical approximation (non-extended cr- 
model) we may re- write the constraint Eq. © in the form: 



Q 2 (r;r';p) = 7 



(16) 



for any value of arguments. The transformation from 
Eqs.© and © to Eas. lfbljl and (fl~5f) is exact, while the 
form Eq. Ijlfijl of the original constraint Eq. © is appropri- 
ate only for the sub-manifold £q. On this sub-manifold, 
the saddle-point equation l|13|) looks like 



V r [Q(r,r';p)V r Q(r,r';p)] = 
sC p /(nPD)[A z k, Q(r', r'; p)]5(r - r') 



(17) 



Similar to [6j , we take into account only the non-compact 
boson degree of freedom (the angle 9; < 8 < oo) and 
parametrize the boson sector of Q in the form: Q RR = 
-Q AA = cosh#, Q AR = -Q RA = sinhfl. In this approx- 
imation, the distribution function V(y\ determined from 
Eq.© by the inverse transformation, reads: 

V{v) = J ds/{2ni)d6exp(sSv/D - F b [6]) , (18) 

where the s-integration goes along a contour parallel to 
the imaginary axis and the boson action F b is 



- ttvD ^— \ f dvdr'„ „, 

*w] = — £ / -T^ [v ^ (r;r;p 



2 ^ I L a 
p 



(19) 



/dr 
-^[cosh0(r;r;p)-l]. 
p 

Integration over s in Ea. (|18f) leads to the constraint: 

1 = 7^^/ ^[cosh0(r;r;p)-l]. (20) 
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The saddle point equation for the action F& Eq . (|19l is 

V 2 r) 0(r; r'; p) - k p sinh%'; r'; p)$(r - r') = , (21) 

where n p = 2sC p /(tti'D). This equation differs substan- 
tially from that derived in 6f for the usual cr-model. 
Ea. (|2*T|l is easy to solve. We find: 



0(r;r';p) = -n p V{r; r') sinh6»(r'; r'; p) . 



(22) 



where — 2?(r;r') is the Green's function of the Laplace 
operator with the boundary conditions: D = on con- 
tacts and V n 2? = elsewhere. The value of 0(r;r;p) is 
determined by the self-consistency condition for Eq. l|22|) 
at r = r' 0] . One of solutions to this equation is trivial: 
0(r;r;p) = 0. Assuming n p <C 1, we find with the log- 
arithmic accuracy a non-trivial solution: 0sign[Re(#)] w 
ln[— 1/{k p T>{t; r))]. To minimize the saddle-point action, 
we take into account the non-trivial solution for only one 
of possible p-channels (and put |p| — pf) - this is sim- 
ilar to the restriction to only one "pre-localized state" 
in |(|. Finding the quantity s from Ea. (|2L)t . we obtain 



0(r;r;|p| = pp) — \n[5v / (vt/S)]. Using Eqs. 
(113, we get finally for V{v), Eq.®: 

V{v) oc exp (-[irvD/{2V)\ In 2 [5v/{vtA) 



and 



(23) 



where 1/2? = L~ d J dr/2?(r; r). The "anomalous" 
saddle-point contribution Eq.J33J) to Viv) has been de- 
rived under the condition 8v/v 3> rA (large 0(r;r)). 
This contribution dominates over the ordinary Gaus- 
sian background V(y) oc exp [— {5v) 2 /{2 < {5v) 2 >)] at 
8v/v > g-^il/L)^ 2 - 1 (we define g as g = 2nE c /A). 
For d = 2, V ~ (27T)" 1 \n(L/l) and we have for Eq. Q : 



2-'(i / ) oc exp 



*9 



21n(L/Z) 



DtA 



(24) 



Eq. (|24|l coincides with Eq. Q (up to the factor 2 resulting 
from the difference between the orthogonal and unitary 
ensembles) and with the result of the physical reasoning 
0, El based on Eq.©, but it has weaker restrictions 
on its validity range. However, at d=3 there is no agree- 
ment between Ea. (|23|l and the result based on a 
version of Eq. (JSJ for d=3 space . Also, there is some 
discrepancy (although not too dramatic) in the quasi- 
one-dimensional case: in [ll|, the logarithm's argument 
contains an extra factor grA ~ E c t. A reason of the 
discrepancies at d ^ 2 is not clear yet. It may be a sign 
of importance of fluctuations and interaction of boson 
and fermion modes, ignored in the simplified saddle-point 
analysis of |6( and in the present paper (fluctuations can 
also restrict the possibility of dealing with only £ sub- 
manifold of the saddle-point manifold £). Note that the 
RG results, Eqs.10 and (jHJ, refer formally to only d=2. 

To conclude, an exact functional integral representa- 
tion, Eq.@, is suggested for the global electron DOS dis- 
tribution function V(y). This representation allows one 



to perform the disorder averaging and to derive a "slow" 
functional, Eqs. 10 and Q. The formalism reproduces 
the diagrammatic results for moments of DOS and pro- 
vides also a possibility to study non-perturbative regimes. 
In particular, the logarithmically normal asymptotics of 
V{v) has been derived. With some complication of the 
functional integrals, the approach can be implemented 
also for the conductance distribution function V(g) 
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